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Abstract 

We analyze the nucleon axial-vector coupling to two loops in chiral perturbation theory. We show 
that chiral extrapolations based on this representation require lattice data with pion masses below 
300 MeV. 
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1. The axial-vector coupling constant is a fundamental property of the nucleon that can e.g. be 
determined in neutron /3-decay (for a review on the nucleons axial properties, see [1]). It is directly 
related to the fundamental pion-nucleon coupling constant by the Goldberger-Treiman relation and 
thus of great importance for the problem of nuclear binding. In the last few years, first attempts to 
calculate qa using various approximations to QCD on a discrctizcd space-time (lattice QCD) have 
been published, see e.g. [2-5]. These results are obtained for quark masses considerably larger than 
the physical ones, the lowest quark masses considered e.g. in the most recent study [5] correspond to 
a pion mass of about 350 MeV (which should already be close to the so-called chiral regime). It is 
therefore necessary to perform a chiral extrapolation to connect these lattice results with the physical 
values of the quark masses.*"^ Long before the advent of these lattice data it was noted that the 
chiral expansion of the axial-vector coupling does not show the expected convergence behaviour for 
an SU(2) quantity — the correction of order is of the order of 30% at the physical pion mass 
although it is two orders down compared to the leading term [8]. One therefore can not expect the 
one— loop representation to be very accurate for increasing pion mass. In fact, the complete (fourth 
order) one-loop result is steeply rising with growing M^r while the lattice data show essentially no pion 
mass dependence A possible solution to this problem was offered in Ref. [9] where an effective 
field theory with explicit delta degrees of freedom at leading one-loop order could lead to a flat pion 
mass dependence of gA, requiring, however, a fine tuning of certain low-energy constants. For a recent 
update, see [10]. Two remarks on that result are in order: First, it should also be noted that most 
lattice data available at that time arc far outside the regime of applicability of the effective field 
theory. Second, to judge upon the usefulness of such an approach requires a systematic analysis of 
many other observables which has not been done so far*^ (for further discussions on this issue, see 
e.g. [11,12]). To obtain a deeper understanding of the chiral expansion of gA it appears therefore timely 
to go beyond the one-loop approximation in nucleon chiral perturbation theory. This is precisely the 
issue of this letter. Using renormalizations group methods, we will determine the coefficient of the 
double log in gA - that arises first at two-loop order - from the existing one-loop results (note that 
the spectral function of the axial form factor to two loops was already worked out in Ref. [13]). We 
also give the general structure of the two loop representation of gA, and determine the numerically 
leading contributions to the single logarithm and polynomial terms at order and M^, generated by 
graphs with insertion proportional to the large dimension two low-energy constants C3 and C4 and the 
dimension three LEG die. We thus achieve more detailed information on the quark mass expansion 
of the nucleon axial-vector coupling constant. 

2. Our calculation is based on the effective Lagrangian of pions and nucleons coupled to external 

sources. The various contributions to S-matrix elements and transition currents are organized in 
powers of the small parameter q, where q collectively denotes small pion four-momenta, the pion 
mass and baryon three-momenta. The effective SU(2) Lagrangian is given as a string of terms with 
increasing chiral dimension, 

where the ellipsis denotes terms not needed in what follows. The local operators at the various 
orders are accompanied by low-energy constants (LECs), these are denoted as q, dj, Cj, . . . for the 

*''ln addition, one has to correct for finite volume and finite size effects, which we do not consider in what follows. For 
recent studies, see [6,7]. 

*'One should be somewhat cautious to draw too strong conclusions from such observations because most of the available 
lattice QCD results are far outside the range of applicability of chiral perturbation theory or any model-independent 
scheme. 

*^This was to some extent attempted in [9] where one LEG was constrained by matching to the pion-nucleon theory, 
in which this LEG had been determined earlier. 
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dimension two, three, four, . . . pion-nucleon Lagrangian and Zj for the masonic LECs of dimension 
four. According to the power counting, the tree approximation is given by tree graphs with insertions 
from The one-loop approximation contains further tree graphs with insertions from C^j^^ and 

one-loop graphs with insertions from C^^j^ and at most one insertion from C^li- At two-loop order, we 
have two-loop graphs with insertions from C^^j^ and at most one insertion from C^}^, one-loop graphs 

with insertions from C^^^ and further tree graphs related to C^j^^ (and the corresponding mesonic 
contributions). Since we are interested in the quark mass expansion of the axial- vector coupling g^, 
it is most convenient to work in the heavy baryon framework (for a review, see [14]). In two-flavor 
chiral perturbation theory, the quark mass expansion is mapped onto an expansion in the pion mass, 
whose physical value is denoted by Mt^. Consequently, the chiral expansion of takes the form 

3A = so |l + (7j^ln^+fti) Mj + QsMj 

|l + A(2) + a(=^) + A(^) + A(5)| + O(M^) , (2) 

with go the chiral limit value of gA, dA = S'oll + C'(M^)], A is the scale of dimensional regularization, 
and A^"-* denotes the relative correction at order M". Further, F denotes the pion decay constant 
in the chiral limit, = F[l + 0{M^)]. To the order we are working, we require the quark mass 
expansion of F^r, 



90 



F^ = F 

in terms of the scale-independent LEG 



/Vf2 



(3) 



167r^£^(A) -21n(M^/A) , (4) 



where ^4 (A) is the corresponding scale-dependent renormalized LEG. Note that this explicit quark 
mass dependence of Ft^ has to be accounted for when one studies the axial coupling as a function of 
the pion mass. Note also that when we generate the numerical value of ^4(A) from ^4, we have of 
course to use the physical value of the pion mass in Eq. Q. Furthermore, the chiral expansion of the 
pion decay constant generates contributions to 04, /34, 74 and 05. This can be seen from Eq. which 
is expressed in terms of the chiral limit value F instead of the physical value, as it is commonly done. 
At a fixed pion mass, these two representations are of course equivalent. The explicit expressions of 
these additional quark mass dependent terms are given below. The third order one-loop coefficients 
a2 and (32 in Eq. ^ were first given in [15] and the one-loop fourth order calculation was completed 
in [8] with (we use the by now standard notation of Refs. [14, 16]) 

as = -2 - 45^ , 
132 = -(d{e{X)-2god^2sW^ 



"3 

with 



goV'"' --—-V (47rF)2 ' 

^ (3 + 3^0 - 4moC3 + SniQC^^ , (5) 



247rF2mo 



rl^ (\\ - J I go(4-go) . 
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As in [16], we set d2s = in what follows. Again, note that these relations for d[(A) {i = 16,28) 
have to be taken at the physical value of when it comes to pin down their numerical values. The 
dimension two LECs 03,04 can be determined e.g. from the analysis of elastic pion-nucleon scattering 
and the dimension three LEG die from the reaction vrA^ — s- ttttN (for a detailed discussion see e.g. 
Ref. [11] and references therein). In Eq. (jSJ, niQ denotes the nucleon mass in the chiral limit. To the 
order we are working, it is related to the physical nucleon mass mi\f via 

TUN = mo- AciMl + 0{Ml) , (7) 

with ci another dimension two LEG that can be determined e.g. from low energy pion-nucleon 
scattering data or the pion-nucleon sigma term. This quark mass dependence of the nucleon mass 
induces corrections at fifth order from the third order coefficient 03 ~ I/ttiat. In what follows, we 
always work with mg and absorb this induced contribution in the combination of LEGs contributing 
to 05. 

In this letter, we are going to evaluate the coefficient 04 of the double logarithm which arises at 
two-loop order. This requires only parameters from the one-loop calculation, as already stressed by 
Weinberg in his seminal paper [17]. The coefficients and as contain combinations of LEGs 

from C^^^"^ and unknown LEGs from C^^'^\ we will estimate these using naturalness arguments 
and also from the description of the available lattice data (at small enough pion masses). Note that 
we can, in addition, work out the numerically large contributions to these coefficients ~ Cj/mo and 
~ Cj/mg from the expansion of the corresponding relativistic one-loop graphs (together with the 
induced contributions from the quark mass expansion of F.,^). 

3. The application of renormalization group (RG) methods to chiral effective Lagrangians was pio- 
neered by Weinberg [17]. He showed that the coefficient of the double log ~ In^ can be entirely 
expressed in terms of coupling constants of the one-loop generating functional. For recent applications 
of such RG methods in chiral perturbation theory for mesons, see e.g. [18-20], a nice discussion of 
this and related issues is given in [21]. Here, we wish to apply the same arguments to the effective 
pion-nucleon Lagrangian. According to the power counting, the double logs are generated from two- 
loop graphs at 0{q^). Employing a mass-independent renormalization scheme (here: dimensional 
regularization) , the two-loop divergences take the generic form 



(8) 



with d the number of space-time dimensions, e = d — A and k{d) is a function of d that depends 
on the specific diagram under consideration. This function can also be expanded around d = 4, 
k{d) = ko + kie + 0{e^). The leading term in this expansion generates the non-local divergence 
~ ko In M/ € that must be cancelled by one-loop graphs with insertions from the dimension three 
effective pion-nucleon Lagrangian (parameterized by the unrenormalized LEGs di). Such graphs give 
the generic contribution 



hi{d) 



Ki , K,^^M^ ^ (4vr)2d[(A) , ,,„,2^rMM. 

A 



2 + - In ^ + ^^-^ + (47r)^ <(A) In ^ + 



2 (47r)4 

where the di are the dimension three LEGs that have the form [16] 



di{d) = X' 



(47r)2e 



+ <(A) + ... 



(9) 



(10) 



where we use the basis of operators enumerated in [16] with the corresponding /^-functions listed 
there. Here, hi{d) is a function specific for the coefficient under consideration, that itself depends on 
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Figure 1: Topologies of the one-loop graphs that generate the coefficient of the double log at two-loop 
order. The hatched square denotes a dimension three insertion proportional to some of the LECs di. 



d via hi{d) = /ijo + hue + 0{e^). As shown by Weinberg, the elimination of the non-local divergence 
is guaranteed by the RG condition 

fco = ^ hiO Ki . (11) 

Using this equation, we can now calculate the coefficient of the double log. There are two types of 
diagrams contributing, namely irreducible and reducible ones, the latter being related to wave function 
renormalization. It is important to note that the notion of reducibility here refers to the two-loop 
graphs. The non-vanishing irreducible two-loop contribution is generated from the graphs shown in 
Fig. ^ The following operators (given in terms of their LECs) contribute to the various graphs: a) 
di6, c?25 b) dio, dii, di2, c?i3, and die, c), di, ^2, du, ^26 and dso d) die and ^29, e) ^24 and ^28, f) 
di6; and ^29, g) ^261 ^27 and ^281 and h) ^24 and d2s- Using the /3- functions from [16] gives the 
contribution to the double log generated by these diagrams. It reads 

af = 2 (1 + - gi) . (12) 

Furthermore, there are reducible graphs generated from wave function renormalization. These are 
given by g]^^°°^ ■ 2'^~^°°p and their contribution to the double log is 

a?^ = 9(5o' + 2ffo') . (13) 

Putting pieces together (i.e. the contributions from the irreducible diagrams, Eq. (|12() . the reducible 
graphs, Eq. (|13j) . and the induced term a^, see Eq. p7|) below), we have thus for the coefficient fco = 04 
of the double log 

a, = aT + aT" + 04 = -y - y <7o' + IG^o' ■ (14) 

This is the central result of this paper and allows us to analyze the leading two-loop correction to 
the axial-vector coupling constant. In the formulation using the Lagrangian given in [16], one has to 
deal with a large number of equation of motion terms. These can be, however, eliminated from the 
effective Lagrangian as done in [22]. We have therefore also performed the calculation in the basis 
given in that paper and employing the pertinent /3-functions. We find the same result as in Eq. (|14|l . 
which serves as an excellent check on our calculation*"^. Note that this procedure generates also part 
of the single log coefficient 74 in Eq. We only give the contribution generated from the operator 
proportional to the axial LEG die (note that some of the other LEGs di are also known, see Ref. [24], 
but only die plays a prominent role in the chiral expansion of g^) 

74'^« = -12die(A)(^ + 5o) • (15) 

*^Provided one corrects for the typographical error in /3ii in that paper, see also [23]. 
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Note further that using a relativistic formulation (see e.g. Ref. [25]), it is easy to work out the 1/mo 
and l/rriQ corrections to the large contribution 03 ~ 94 (for the parameters given below, see also the 
discussion in [11]) that will give a sizeable contribution to the coefficients /34,74 and as, respectively. 
These terms are given by 



74 



4(C4 - C3) 

mo 



mo 47r^F-^ ' 

= ;|i6^- ^''^ 

The numerical values of these contributions will be given below, but we remark already that in 
particular 74' will contribute sizeably. Finally, we collect here the induced terms form the quark mass 
expansion of the pion decay constant, cf. Eq. Q. Denoting these by a tilde, they read 

d4 = 4a2 , 

2 4^2 
74 = --p^a2l4W 



- (47rF)V2^^^^) 



as 



We end this section with a brief comment of the one-loop chiral EFT representation given in [9]. 
Including an explicit delta to leading one-loop order generates some of the terms ~ and part of 
the coefficients 74 and P4 (plus some other higher order terms ~ M^^/A^™ with 2n — 2m = 2 and A is 
the delta-nucleon mass splitting). However, already at third order in the pion mass, this neglects other 
resonance contributions to the LECs C3 and C4 (for a detailed discussion, see e.g. [26]) and therefore 
that representation can not be considered as accurate as the one developed here (for pion masses in 
the chiral regime). 

4. We are now in the position to put pieces together. First, we consider the contribution of the various 
(incomplete) terms at fourth and fifth order in the pion mass for the physical values of the quark 
masses. For that, we express all parameters in terms of their chiral limit values. We take £4 = 4.33 
corresponding to F = 87MeV and mo = 880 MeV. If not stated otherwise, we use diQ = —1.76GeV~^, 
C3 = 3.5 GeV~^ and C4 = — 4.7GeV~^. We also work at A = mo- Note that we have varied these LECs 
within their allowed ranges, but this did not lead to any sizeable changes to the results given below. 
Furthermore, we vary go between 1.0 and 1.2. To be definite, let us set go = 1. Collecting pieces, we 
obtain 

= af +'^"'^ + d4 = 31 - 24 = 7 , 
74 = + 74 + 74^« = (37.3 + 3.2 + 74.8) GeV-^ = 115.3 GeV"^ , 
/34 = + f34 = (13.3 + 0.9) GeV"^ = 14.2 GeV"^ , 

as = ag' + ds = (-16.0 - 4.3) GeV"^ = -20.3 GeV"^ . (18) 

Of course, the coefficients 74, and as receive further corrections from LECs that have to be de- 
termined e.g. from an analysis of lattice data or estimated assuming naturalness. We denote these 
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Figure 2: The axial-vector coupling as a function of the pion mass. Solid (red) line: qq = 1.2, dig = 
-1.76GeV-2, -fl = 50GeV-2, (3( = eOGeV-^, a{ = 20GeV-5. Dot-dashed (black) line: go = 
1.1, di6 = -0.92GeV-2, 7/ = 40GeV-2, /?{ = 20GeV-^, a{ = 50GeV~5; Dashed (green) line: 
go = 1.0, die = -1.76GeV-^ 7/ = -SOGeV'^ = = 0. The dotted (violet) line is the complete 
one-loop result with go = 1, diQ = —1.76 GeV~^ and using the physical values of the nucleon mass and 
the pion decay constant. The (magenta) circle denotes the physical value of gA at the physical pion 
mass and the triangles are the lowest mass data from Ref. [5]. 

additional contributions by 7/ , /jf and al- For the moment, we set 74 = /^f = = 0. In the notation 
of Eq. ((2)), these results translate into 

A(2) = -15.3 % , 

A(3) = 25.6 % , 

^(4) ^ ^(4) ^ ^{4) ^ ^(4) ^ ^Q^g _ ^ Q ^ _5^g ^ 

^(5) ^ % _ (19) 

The sizeable fourth order contribution is entirely due to the large coefficient 74, largely due to the 
insertion of the operator ~ dig, cf. Eq. (|15j) . Note also that the fifth order term is very small at the 
physical point. Ignoring the higher order corrections, one can calculate the chiral limit value of gA 
from Eq. using the values collected in Eq. (|19|) . 

go = 1.21 [1.12] , (20) 

where the number in the brackets refers to the choice Jig = — 0.92GeV~^ and we use gA = 1.267. 
Of course, these numbers will be affected by the unknown LEG contributions 7/, /Sj and a^. These 
values are consistent with the findings in [7] . We also remark that the chiral expansion oirriN is much 
better behaved and one thus can successfully apply one-loop extrapolation functions to pion masses 
below 450 MeV (for detailed discussions, see e.g. [11,27,28]). 

We show in Fig. 121 some typical examples for the pion mass dependence of gA for values of 7/ , /jf , al^ 
that lead to an approximately flat behaviour for not too high pion masses. These values are of natural 
size as a comparison with the induced pieces collected in Eq. (|18j) reveals. This is very different from 
the one-loop representation, which fails to generate a flat quark mass dependence for values above 
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the physical pion mass, see e.g. [9,11]. We have studied many more combinations of the LECs and 
found that for this representation to be useful (that is leading to a moderate theoretical uncertainty), 
the pion mass should be less then 300 MeV. This can also been seen if one compares to the complete 
one-loop result as depicted by the dotted line in Fig. I2j The existing lattice results are at still too high 
pion masses for a model-independent extrapolation to the physical values of the quark masses. For a 
particular choice of the LECs, we can describe the trend of the lattice data up to Mj^ ~ 600 MeV, but 
the theoretical uncertainty is simply too large for such values of the pion mass. 

5. In this letter, we have studied the pion mass dependence of the nucleon axial-vector coupling 
constant g^. This is a fundamental observable for our understanding of the nucleon structure in 
the regime of strong QCD. First lattice simulations have appeared and so far, chiral extrapolation 
functions appearing in the literature are either based on (leading) one-loop chiral effective field theory 
results with explicit deltas (for a critical discussion, see e.g. [12]) or are very model-dependent. We have 
provided the two-loop representation in baryon chiral perturbation theory, see Eq. (jSJ , and determined 
the coefficient of the double log term ~ M^ln^iW^ based on renormalization group arguments. We 
have also determined some numerically important contributions to the terms ~ In M^^, and M^. 
We have shown that with LECs of natural size one can indeed obtain a flat pion mass dependence of 
gA for pion masses below 400 MeV. We conclude that lattice data for pion masses below 300 MeV 
are required to use this representation with a moderate theoretical uncertainty. Such data should be 
available in the near future. 
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